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Abstract 

We consider the influence of realistic island diffusion rates to homoepitaxial 
growth on metallic surfaces using a recently developed rate equation model which 
describes growth in the submonolayer regime with hyperthermal deposition. To this 
end, we incorporate realistic size and temperature-dependent island diffusion coef- 
ficients for the case of homoepitaxial growth on Cu(100) and Cu(lll) surfaces. We 
demonstrate that the generic features of growth remain unaffected by the details 
of island diffusion, thus validating the generic scenario of high density of small is- 
lands found experimentally and theoretically for large detachment rates. However, 
the details of the morphological transition and scaling of the mean island size are 
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strongly influenced by the size dependence of island diffusion. This is reflected in 
the scaling exponent of the mean island size, which depends on both temperature 
and the surface geometry. 
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copper 



1 Introduction 

Hyperthermal deposition (HTD) techniques, such as ion-beam assisted depo- 
sition (IBAD) and low-energy ion deposition (LEID) [1] have recently been 
shown to have great potential in controlling and improving the properties of 
thin films as grown by molecular beam epitaxy. In HTD the island density 
is larger, the average island size is smaller [2,3], and island size distributions 
are considerably broader [1,3] than in ordinary thermal deposition. Possible 
atomistic processes responsible for these effects include ion enhanced mobili- 
ties, cluster dissociation [1], and defect creation during deposition [2,4]. 

A particularly striking observation made in the LEID experiments is that with 
different deposition energies an anomalously high density of small islands is 
observed, with the scaled distribution function behaving as f(x) ~ 1/x for 
x < 1 [1]. We have recently shown using the rate equation (RE) approach 
[5,6] that this anomalously high density of small islands is due to a unique 
balance between island-island aggregation and enhanced adatom detachment. 
These studies were aimed at describing the generic features of HTD, and thus 
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relatively idealized approximations for the terms in the rate equations were 
used to obtain analytical estimates for the relevant growth exponents and 
the scaling function. One of the most important questions that still remains 
open is the role of island diffusion, since for mobile islands the aggregation 
rates in the RE approach depend explicitly on the diffusion coefficients D s 
of islands of different sizes s. Island diffusion on surfaces has been studied 
both theoretically [7,8,9] and experimentally [10,11,12,13,14,15,16]. While in 
the large island limit the size dependence of the island diffusion coefficients 
can be classified based on simple basic processes [17], for smaller islands D s 
depends on the geometric and energetic details of the underlying surface, and 
can be a complicated, non-monotonic function of s [18,19]. 



In this work, our aim is to explore in detail the influence of realistic island 
diffusion coefficients to submonolayer growth with HTD. To this end, we em- 
ploy the RE model of Ref. [5] and replace the usually assumed idealized power 
law forms of D s with realistic, temperature and size-dependent diffusion co- 
efficients D a (T) for Cu on Cu(100) and Cu(lll) surfaces. These two systems 
highlight the large differences which occur for surfaces with different geome- 
try and energetics. We demonstrate that while the scaling function of the size 
distribution is largely unaffected by the details of the diffusion coefficients, the 
quantitative values of the growth exponents are sensitive to island diffusion. 
These predictions can be easily tested by HTD experiments on different Cu 
surfaces. 
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2 Model 



In HTD the reversibility of growth is manifested through enhanced adatom 
detachment from islands [1,4]. Detachment and island mobility allow us to ne- 
glect spatial correlations between islands [20,21,22,23], which justifies the rate 
equation description of the problem. Thus, growth is driven by the interplay 
between aggregation and detachment, and can be schematically expressed to 
be composed of reversible events Ai + Aj — > A i+ f, Aj — > A\ + A,-_i between 
islands of sizes i and j with the rates of aggregation and detachment specified 
by reaction rates K(i,j) and F(i,j), respectively. The corresponding REs for 
the areal density n s of islands of size s > 1 read as [23,24] 

-JT = \ J2 [K(h3)mnj - F(i,j)n s ] 

m Z i+j=s 

oo 

- Y\ K ( s ij) n s n 3 ~ F(s,j)n s+j ] + $6 la , (1) 

where $ is the deposition flux of adatoms in monolayers per second (ML/s). 

The aggregation rate K(i,j) for islands of sizes i and j with diffusivities D i 
and Dj is given by the Smoluchowski formula K(i,j) oc (Di+Dj), which is also 
consistent with the point island approximation used here [20,24]. Previously 
[5] we used a power-law form for the diffusion coefficients Di ~ with 
1 < fi < 2 appropriate for island diffusion on metal surfaces. In order to 
study the effects of details of the aggregation rates in observable quantities 
we replace the idealized power-law form by realistic size and temperature 
dependent diffusion coefficients, as discussed below in more detail. 

The detachment rate of adatoms from islands of size % + j = s depends on 
the island size, but only detachment of single adatoms is allowed. Thus the 
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detachment rate is given by F(i,j) = F (i + j) a (Su + Sij) , where the exponent 
a is in principle a parameter, but is chosen to be a = 1/2 in the present study. 
In LEID this form for the detachment rate is physically a well justified choice, 
because in the regime of bombarding energies from 10 eV to 100 eV adatom 
detachment dominates and island breakup into larger pieces is not expected 
to occur [4,25]. Moreover, since every deposition event at the vicinity of an 
island boundary can be assumed to detach adatoms at least with a probability 
proportional to the island perimeter (i.e. s 1 ^ 2 ), a = 1/2 is a reasonable lower 
limit. This has been also confirmed by recent Molecular Dynamics simulations 
on ion bombardment enhanced detachment in island size region up to 25 atoms 
where values of 0.4 < a < 0.6 were found [25]. 

We solve the rate equations using the particle coalescence method (PCM) 
[22,20,24]. PCM employs a point-island approximation, which is valid at low 
coverages or large island separations. In PCM aggregation and detachment 
events occur with probabilities specified by the corresponding reaction ker- 
nels, and the deposition with the rate proportional to the given adatom flux. 
An event is then randomly chosen with a probability weighted by the corre- 
sponding rate. Since REs describe the system in the mean-field limit, there is 
no information on spatial correlations in the system. Therefore, in the simula- 
tions it is sufficient to construct a list of islands, which does not correspond to 
a physical lattice. To conduct an aggregation event, for example, two islands 
are randomly chosen from the list, and an attempt to aggregate them is made. 
The complete mixing of the islands required by the mean-field approximation 
[22] is thus implemented much faster than including island jumps into empty 
lattice sites [24]. 
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3 Diffusion Coefficients 



In reality the diffusion coefficients of islands depend on island size non-trivially, 
and not simply as power laws. They can even oscillate as a function of the 
island size [10,14,8]. These oscillations can be in part understood by the sur- 
face geometry: certain close packed island configurations can be more stable 
than others on certain surfaces [19]. Here we concentrate on two qualitatively 
and quantitatively different types of diffusion coefficients, namely those of Cu 
adatom islands on Cu(100) and Cu(lll) surfaces. 

First, the diffusion coefficients D S {T) from kinetic Monte Carlo (KMC) simula- 
tions, which were based on effective medium theory energetics for the Cu(100) 
surface, are shown in Fig. 1 at three different temperatures T = 300, 500, 
and 700 K [8]. At the lowest temperatures the oscillations as a function of 
s are clearly seen. As temperature increases, the oscillations dampen, and at 
the highest temperature the diffusion coefficient is rather well described by 
-D(s) ~ s 3 / 2 , if island diffusion is limited by atomic motion along the perime- 
ter of an island [17]. For Cu(100), there is a strong dependence of D s on 
temparature here because the single adatom diffusion barrier is about 0.4 eV. 

In contrast to the behavior on Cu(100), island diffusion on Cu(lll) follows 
a completely different trend. The most spectacular difference is the almost 
complete absence of oscillations in the diffusion coefficient with size for the 
smaller islands (1 — 10 atoms) [26]. This is seen to occur when small islands 
diffuse via occupancy of both hep and fee sites. Also, there is relatively little 
temperature dependence on D S {T) for small sizes, since the surface is rela- 
tively smooth with an adatom activation energy of 0.026 eV, and almost the 
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Fig. 1. Size and temperature dependence of diffusion coefficients for 2D adatom 
islands on Cu(fOO) at three different temperatures from the KMC simulations of 
Ref. [8]. For the largest island sizes the data approximately follow the power law 
D s ~ s 3 / 2 . 

same for dimers and trimers. In the case of Cu(lll), the calculations were 
carried out using a KMC procedure in which the system is allowed to evolve 
according to diffusion processes of its choice. This is facilitated through the 
automatic generation of a database of possible diffusion processes and their 
activation energy barriers calculated using embedded atom method poten- 
tials. The storage and retrieval of information from the database is done via 
a pattern recognition scheme [27]. From simulations performed at T = 300 K, 
500 K, and 700 K, the effective diffusion barriers are found to increase almost 
monotonically with size, while the diffusion coefficient takes the form shown in 
Fig. 2. Calculations performed for larger sized islands (19 — 100) show a power 
law scaling of D with island size and with an exponent of about 1.57 which 
is mildly temperature dependent [28]. This is consistent with the theoretical 
value of 3/2 for large island motion dominated by atom diffusion along island 
perimeter [17]. 
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Fig. 2. Size and temperature dependence of diffusion coefficients for 2D adatom 
islands on Cu(lll) at three different temperatures from the KMC simulations. The 
dashed lines indicate fit to a s~ 3 / 2 power law [17]. 

4 Scaling 

In our previous study [5] the island size distributions were of the scaling form 
and the mean island size had a power-law form. Therefore, even if the de- 
tachment rate was not a homogeneous function of the island size, we could 
extract well-defined effective scaling exponents for the mean island size and 
the size distribution function. In the present case, we expect that scaling is 
invalidated due to realistic diffusion coefficients which have a naturally inho- 
mogeneous form. There is no guarantee of either the existence of usual scaling 
type of solutions or uniquely defined scaling exponents for the mean island size 
and island density. Thus, e.g. the dynamic exponent for the mean island size 
can be expected to be model-dependent, i.e. we could obtain different values 
for different sets of diffusion coefficients. This has been shown to be the case 
in more general reversible growth models [29] and growth with mobile islands 
[20,30], where the scaling exponents explicitly depend on the details of the 
model. 
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We define the mean island size s = M 2 /Mi, where the k th moment of the 
distribution is defined as = J2 s >i s k n s . If the scaling of the island size 
distribution holds, we expect that n s = 9s~ 2 f(s/s), where s ~ 9 13 , 9 is the 
surface coverage, and f(x) is a scaling function independent of coverage. In 
our previous work [5] we found that the scaling function is always singular, 
i.e. f(x) ~ 1/x for x < 1. In order to compare the scaling function with our 
previous results, we will use in the following the modified scaling function of 
the form g(x) = xf(x). 

5 Results 

The PCM simulations were carried out using island diffusion coefficients dis- 
cussed above for Cu(100) and Cu(lll), and with the detachment rate char- 
acterized by the exponent a = 1/2. We define k = F /Di and 1Z = -Di/<E>, 
where D 1 is the adatom diffusion coefficient. Parameters k and R denote the 
importance of detachment relative to diffusion, and of diffusion relative to de- 
position flux, respectively. In the simulations the corresponding values were in 
the ranges 1CT 6 < k < 1CT 1 and 10 5 < R < 10 9 . The simulations show that 
for large detachment rates, s and iV indeed follow a power-law type of behav- 
ior, but as k decreases, (3 becomes coverage-dependent, and saturates only for 
k — > as we have previously shown [5]. Thus, only for large k can well-defined 
scaling exponents be extracted and regular island growth observed, and in this 
regime the island size distributions are of scaling form. 

The measured values of the dynamic growth exponent (5 at different tempera- 
tures and large k both for Cu(100) and Cu(lll) are shown in Table 1. The ex- 
ponent (5 seems to be temperature dependent for Cu(100) but not for Cu(lll). 



9 



For Cu(100) this can be explained by examining the corresponding curves for 
the diffusion coefficients in Fig. I. For small sizes (which is the case for large re) 
the average slope of the diffusion coefficient curves depends strongly on tem- 
perature. If one fits a power law into the initial part of the data (for s < 20), 
the effective exponents are /i e s ~ 5.4, 3.6, and 2.8 at T = 300, 500, and 700 K, 
respectively. For the power-law type aggregation kernels we found [5] that the 
dynamic exponent behaves as (5 — l/(a + (i), when K(i,j) oc Using 
the effective exponents above this prediction gives (5 ~ 0.17, 0.24, and 0.30, 
showing a similar trend as the measured values. If the measured values are 
used, somewhat smaller /x e fr are obtained than from the fits to the diffusion 
coefficient data. If the fitting is done only through maxima of the diffusion 
coefficient curves, a better agreement is obtained. It is also interesting to note 
that s does not depend on 1Z in this regime, but only on re. This suggests that 
it could be possible to tune the regime for the effective diffusion exponent on 
Cu(100) by changing re. However, this effect is probably too small to be seen, 
since e.g. on Cu(100) at T = 300 K we get s ~ 12 atoms for re = 10~ 3 , and 
s ~ 7 atoms for re = 10~ 2 , while the differences in j3 between these two cases 
are within the errors. For Cu(lll) there is no clear power-law for small sizes 
where the simulation data for diffusion coefficients exist. Instead, small island 
sizes seem to be rather mobile relative to adatoms in all temperatures. 

We expect to observe scaling of island size distributions on the basis of the 
fact that well-defined scaling exponents exist for large detachment rates. In 
Fig. 3 we show scaled island size distributions for large values of detachment 
rates using the diffusion coefficients for Cu(100) and Cu(lll) shown in Figs. 1 
and 2, and setting a = 1/2, K = 10 6 , and re = 10~ 2 . On Cu(100) the data 
for different temperatures collapse to a single curve, and for x < 1 the scaled 
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T (K) 300 500 700 
Cu(100) 0.25 0.34 0.43 
Cu(lll) 0.25 0.22 0.29* 

Table 1 

Dynamic scaling exponents j3 as denned in the text at different temperatures ob- 
tained from the PCM simulations, using the diffusion coefficients for Cu(100) and 
Cu(lll) in Figs. 1 and 2. The asterisk for Cu(lll) at T = 700 K indicates that the 
exponent is not yet constant in time. The errors in other cases are about ±0.05. 




Fig. 3. The scaled island size distributions on Cu(100). The different symbols cor- 
respond to T = 300 K (circles), T = 500 K (squares), and T = 700 K (triangles) 
with TZ = 10 6 , k = 10~ 2 , and 6 < 0.25 ML. The inset shows the distributions on 
Cu(lll) with the same parameters. 

distribution is almost flat. Deviations between different temperatures occur 
at small sizes, which reflects the differences in the diffusion coefficients. The 
inset shows the scaled distribution on Cu(lll). In this case scaling of the 
distribution is not as good, and differences at small sizes are larger than on 
Cu(100). This reflects the fact that on Cu(lll) the diffusion coefficients for 
small island sizes do not follow any power-law type behavior but are almost 
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constant in s, thus invalidating the scaling behavior as in a more simplified 
model for ((J,, a) = (0,1/2) [6]. It is, however, expected that deviations could 
be seen in the large size tail of the distribution. Since the diffusion coefficients 
decrease several orders in magnitude, e.g. from single atoms to islands of size 
10, all aggregation events leading to large sizes are basically between a small 
island and a large one. 

6 Discussion and Conclusions 

To summarize, we have studied a rate equation model with aggregation and 
enhanced adatom detachment corresponding to hyperthermal deposition con- 
ditions using realistic size and temperature dependent diffusion coefficients for 
Cu(100) and Cu(lll). We have shown that qualitatively the generic features 
of growth, e.g. the form of the scaling function of the size distribution, are 
not influenced by diffusion, but quantitative predictions depend on the mi- 
croscopic parameters. For example, the mean island size follows a power-law 
form in all cases with the exponent dependent on temperature and surface 
geometry. Scaling of the size distribution and the mean island size can be 
observed, however, only for large values of adatom detachment rates. These 
findings suggest that differences between various surface geometries can be 
observed only in the dynamic quantities such as the scaling exponent of the 
mean island size, but not e.g. in the scaled island size distribution which is 
time-independent. These predictions can be easily tested with HTD on Cu 
surfaces. 
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